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ABSTRACT 

Several studies discussing imaging polarimetry observations of protoplanetary disks use the so-called radial Stokes parameters and 
to discuss the results. This approach has the advantage of providing a direct measure of the noise in the polarized images under 
the assumption that the polarization is azimuthal only, i.e., perpendicular to the direction towards the illuminating source. However, a 
detailed study of the validity of this assumption is currently missing. We aim to test whether departures from azimuthal polarization 
can naturally be produced by scattering processes in optically thick protoplanetary disks at near infrared wavelengths. We use the 
radiative transfer code MCFOST to create a generic model of a transition disk using different grain size distributions and dust masses. 
From these models we generate synthetic polarized images at 12pm. We find that even for moderate inclinations (e.g., i = 40°), 
multiple scattering alone can produce significant (up to ~ 4.5% of the image, peak-to-peak) non-azimuthal polarization reflected 
in the C/^ images. We also find that different grain populations can naturally produce radial polarization (i.e., negative values in the 
Qtp images). Despite the simplifications of the models, our results suggest that caution is recommended when interpreting polarized 
images by only analyzing the and images. We find that there can be astrophysical signal in the images and negative 
values in the images, which indicate departures from azimuthal polarization. If significant signal is detected in the C/^ images, we 
recommend to check the standard Q and U images to look for departures from azimuthal polarization. On the positive side, signal in 
the U(p images once all instrumental and data-reduction artifacts have been corrected for means that there is more information to be 
extracted regarding the dust population and particle density. 
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1. Introduction 


With the advent of high contrast imaging polarimeters such 
as NaCo/VL T ( [Lenzen et al.||2003]), HiCiao/Subaru (|Tamura| 
et al.||200^ , GPI/Gemini South ( Macintosh et al.||2QQ8| ), and 
SPHERE/VLT ( [Beuzit et al.|2008 l, polarized differential imag¬ 
ing (PDI) has become a standard tool to directly image and an¬ 
alyze the dusty circumstellar environments around young and 
evolved stars (e.g.|Canovas et al. |2012[[2015at [Jeffers et al.|2012 


Garufi et al.|2013[|Min et al.|20i3| [defers et al.|2014| ). In disks, 

PDI have proven powerful to show a variety of features like spi¬ 
ral arms, bright inner rings and/or disk asymmetries revealing 
complex morphologies. However, the yield products of this tech¬ 
nique are generally polarized intensity (Pj) maps, and rarely in¬ 
tensity and fractional polarization (P) maps. The optical proper¬ 
ties of dust grains and their dependence on composition and scat¬ 
tering angle have a significant impact on the polarization of the 
scattered light (e.g., |Murakawa|20IQ[[Min et al.|20I2| ). Multiple 
scattering and optical thickness can also modify the observed po¬ 
larization. Therefore, Pj maps may not always trace accurately 
the underlying surface brightness distribution, and even less so 
the surface density distribution (as in AB Aurigae, see |Perrin| 
let al.|2009D . 


Recently it has become popular to analyze linearly polarized 
images of protoplanetary disks using the so-called radial-Stokes 
formalism ( |Schmid et al.||2006"} |Avenhaus et al.|[MT4l ). This is 
extremely useful because it allows to quantify and control the 


calibration errors potentially to high degrees of accuracy. This 
formalism relies on the assumption that the linear polarization 
in disks is perpendicular to the scattering plane defined by the 
light source (star), the scattering particle, and the observer, i.e., 
it is purely azimuthal. Any non-azimuthal polarization is associ¬ 
ated to noise and/or calibration errors. While this is valid when 
single scattering dominates in nearly face-on disks, the assump¬ 
tion may break down for higher disk inclinations and for peculiar 
dust properties (as for tho se of comets, e.g., jKelley et al.|2004|), 
when grains are aligned ( |Matsumura & Seki||I996[ jWolf et al.| 
|2Q02| ), and/or when disks are optically thick and multiple scat¬ 
tering occurs ( |Bastien & Menard|1988| ). In this letter we explore 
some of these deviations and quantify their impact on the use of 
the radial-Stokes formalism at near-infrared (NIR) wavelengths. 


2. Modeling 

The goal of this study is to demonstrate that departures from 
centro-symmetric azimuthal polarization may occur in proto¬ 
planetary disks for several cases. To prove this we build a sim¬ 
ple but realistic model of a transition disk using th e 3D radia¬ 
tive transfer code MCFOST ( |Pinte et al.|2006[[2009| ). The radial 
structure of the model is similar to that found in the well studied 
disks RXJI 633.9-2442 and Sz91 dCieza et al.|[MT2l |Canovasl 
|et al.|20I5b| ): the innermost region is completely devoid of dust, 
followed by a dust-depleted region extending up to the cavity 
radius (“Region I”), and the outer disk (“Region 2”). 
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2.1. Disk Structure 


The exact choice of the disk geometry and stellar type has a min¬ 
imal impact on our study of the polarization when compared to 
the effect of dust grain properties and/or optical thickness. We 
assume the disk to be at the distance of the Taurus star form¬ 
ing region (140 pc), with central star with mass of IM©, effec¬ 
tive temperature of T* = 4300 K, and radius of = 2.6 Rq. 
This roughl y corresponds to a K 6 star of age ~ 1 Myr ( |Siess 
et al.poQO] ). The disk has an inner cavity of 40 au in radius 
(/?cav) containing 10“^ Mq of dust (optically thin at NIR), with 
the innermost 10 au totally devoid of dust. The surface den¬ 
sity distribution of this inner disk is governed by a power-law 
E(r) = Eioo (^/ 100 au)“^ where Eioo is the surface density at 
r = 100 au. The surface density distribution of the outer disk 
is described by a tapered-edge profile 


E(r) = Ec r ^ exp 


Rc 


2-r 


( 1 ) 


where l^c is the surface density at the characteristic radius Rc, 
and y corresponds to the v iscosity power-law inde x in accre¬ 
tion disk theory (v cx [Hartmann et al.|p~998|). The scale 


height is described assuming that the dust follows a Gaussian 
vertical density profile in both the inner and outer disks H(r) = 
Hioo(r/100mi)^, where ^loo is the scale height at r = 100 au 
and ij/ is the flaring parameter of the disk. The optical depth is 
changed by using different outer disk’s dust masses. 


2.2. Dust properties 

Models show that the polarizability curve of different grain types 
is positive (polarization perpendicular to the scattering plane) in 
most cases, but it can have a negative branch producing polar¬ 
ization parallel to the scattering plane (i..e, radial polarization, 
see [Perrin et al.|[2009[ [Canovas et al.[[2013t [Kirchschlager &[ 
Wolf|2014[ ). Although the detailed composition of disk grains is 
difficult to constrain, observational evidence indicates that cos¬ 
mic dust particles are mostly composed of silicates and, in lower 
quant ity, carbonaceous particles ( [Draine & Lee[p^84l [Draine[ 
2003| ). For our modeling we use grain particles composed of 
a mixture in volume of 70% astronomical silicat es ([Draine & 


Lee[ 1984|) and 30% amorphous carbon particles ( [Li & Gre^ 
berg [1997]) following a standard power-law size distribution 


dn(a) oc a^^-^da using different values for the minimum and max¬ 
imum grain sizes. The full scatteri ng matrix o f the dust popula¬ 
tions is computed with Mie theory ( [Mie|1908| ), using two differ¬ 
ent methods: assuming distributions of 1 ) compact homogeneous 
spheres, and 2) hollow spheres (DHS, [Min et al.[2005[ ) with max¬ 
imum volume fraction /max = 0.8. The minimum and maximum 
grain sizes in Region 1 is fixed and these grains are compact. 
For the outer disk, we explore different values of porosity and 
minimum and maximum grain sizes. The optical indices are de¬ 
rived using effective medium theory (using Bruggeman’s rule). 
The model’s parameters are summarized in Table 


Table 1. Model parameters. Subscripts refer to the two regions of the 
disk (“1” refers to the inner disk, “2” refers to the outer disk). 


Parameter 


Symbol 


Value(s) 


Inner radius 
Cavity radius 
Dust mass 

Surface density index 
Scale height 
Flaring index 
Min. grain size 
Max. grain size 
Porosity 


Rin 

^cav 

^dusti 

P 

^lOOi 

^mini 

^maxi 

P 


10 au 
40 au 

1 X 10“® Mo 
-1 

10 au 
1.15 
3 pm 
20 pm 
0 [%] 


Characteristic radius Rc 

Dust mass Mdust 2 

Surface density index y 

Scale height 

Flaring index \l /2 

Min. grain size 
Max. grain size ^^max 2 


100 au 

0.5 X [10-^ 10“^ 10-4] 
0.3 
10 au 
1.15 

0.03,0.2,5, lOyum 
1000,1500,2000yum 



0 50 100 

Scattering Angle [degs 


50 100 150 

Scattering Angle [degs] 


Fig. 1. a) Polarizability curves computed at 2.2pm using solid spheres 
for amax = 1 mm and different values of (see legend). Solid (dashed) 
lines represent 0% (30%) porous grains. Negative values indicate radial 
polarization, b) Same as a) using DHS theory with /max = 0.8. 


gaussian point spread function (PSF). After convolution, the po¬ 
larized intensity image is constructed as P/ = The 

radial-Stokes parameter^ are defined as 

Qcp = +2 cos( 20 ) -r U sin( 20 ) 

f /0 = -Q sin(20) -r U cos(20) 

where 

0 = arctan--r 6, (4) 

V and y are the coordinates in the image with respect to the star 
(•^ 0 , yo) and 6 is the offset created by instrumental polarization 
(zero in our models). 


( 2 ) 

( 3 ) 


2.3. Synthetic scattered light images 

We use MCFOST to generate synthetic images (Stokes I, Q and 
U maps) at 2.2yum of the disk models for several inclinations 
ranging from pole-on to edge-on. To facilitate the comparison 
with real observations these images are projected into a grid 
with pixel size of 07027 px“^ (equal to the scale of the S27 cam¬ 
era of NaCo/VLT) and convolved with a 2.5-px width (07067) 
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3. Results and discussion 

3.1. The effect of dust properties 

The angle at which incident light is scattered and the polarized 
flux that is produced strongly depends on the detailed properties 

^ The nomenclature/symbols for these quantities is not standardized. 
For consistency with previous studies we use the term “radial-Stokes”. 
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of the scattering dust particles. The polarizability curves of dif¬ 
ferent grain distributions prove this dependency for single scat¬ 
tering of incident unpolarized light. Fig. shows that different 
grain populations can rotate the polarization plane from perpen¬ 
dicular (positive branch) to parallel (negative branch) with re¬ 
spect to the scattering plane. The first case produces azimuthal 
polarization patterns (positive signal in the 2^ images) while 
the second case produces radial polarization (negative values in 
Qcj)). In our models, both the DHS and the solid spheres pro¬ 
duce significant amounts of radial polarization for ^min < 5 pm. 
Using solid spheres, radial polarization is expected only for an¬ 
gles > 130° (backward scattering). For a protoplanetary disk 
this would correspond to scattering on its far side at inclinations 
i > (40° -r a), where a is the flaring angle of the disk. On the 
other hand DHS predicts radial polarization for scattering angles 
<70° (forward scattering) and 140° - 160° (backward scatter¬ 
ing), noting that different values of /max can remove the negative 
branch for forward scattering angles ( |Min et al.|2005| ). 

Radial polarization at large scattering angles (backward scat¬ 
tering) has been ob served at optical and NIR wavelengths in a 
number of comets ( [Kelley et al.||2004[ [Woodward et al.|[2QTT] ), 
and dust grains in protoplanetary disks can have similar com- 
positions to those found in comets ( [Bouwman et al.|2003[ ). Ex¬ 
perimental measurements in the optical show that ensembles of 
arbitrarily shaped dust aggregates als o produce radial polariza- 
tion patterns at large scattering angles ( [Munoz et al.|2000[[2006[ 
[Volten et al.[[2007| ). At optical wavelengths models show that 
porous grains can also produce radial polarization in selected re- 
gions of protoplanetary disks at large inclinations ( [Kirchschlager[ 
& Wolf|2014[ ). Therefore, dust properties alone can produce ra¬ 


dial polarization in particular for backward scattering. 


3.2. The effect of multiple scattering 

Multiple scattering of previously polarized light may lead to 
complex polarization patterns. In particular, the polarization 
plane of the last scattering event could have an orientation that 
is neither parallel nor perpendicular to the scattering plane of the 
first scattering event. This would be reflected in the images. 
To describe this case while minimizing the impact of the grain 
properties, we focus in models having amm > 5pm, i.e., mod¬ 
els with polarizability curves without negative branches (as the 
green curve in Fig. [^. The contribution of multiple scattering 
can be estimated by comparing models with the same grain pop¬ 
ulations but different masses (i.e., different number of scattering 
particles), as higher masses will produce more scattering events. 

Three representative cases are shown in Fig. For each 
model, the images are normalized with respect to the Pj image to 
facilitate a comparison between the and images. All 

models are computed using P = 30%, ajnm,max = 5, lOOOyum, 
i = 40°, and DHS. From top to bottom, the outer disk’s dust 
masses are Mdust 2 = 0-5 x [10“^, 10“^, 10“"^] M©, corresponding 
to integrated optical depths at 2.2 pm along the disk’s mid plane 
of T 2.2 = [10“^, 8.3, 83.0]. By construction there is no instrumen¬ 
tal polarization in our models and the numerical noise from the 
Monte Carlo simulation is at most 0.4% of the Pj image (peak 
to peak, see the optically thin model in top row of Fig|^. 

We find that multiple scattering can naturally produce ob¬ 
servable signatures in the unconvolved U^p images even in disks 
with moderate inclination (i = 40°). This is reflected in the U^p 
images in Fig.|^ which can have maximum values up to ~ 4.5% 
of the Pj image (peak-to-peak). In all our models the signal in the 
U(p image relative to the Pj image increases with inclination and 
with mass. The Q(p images, on the other hand, are mostly insen¬ 


sitive to this effect. Both the DHS and compact spheres methods 
produce the same trends. Additionally we find that increasing the 
porosity slightly increases the signal in f/^, but this feature be¬ 
comes undetectable after convolution with the PSF. Instrumental 
effects like PSF convolution can remove these trends and/or even 
create artificial signatures in the images. 

Models based on Mie theory are an attempt to represent the 
complex, fluffy shapes observed in the dust particles of the inter¬ 
stellar medium. For example, detailed, multi wavelength studies 
of the optically thin debris disks HD 181327 and HR 4796 A ( [Le-[ 
[breton et al.|2Q12t [Milli et al.|2Q15[ ) show that neither by using 
solid spheres nor DHS is it possible to reproduce the observed 
scattering properties of their dust grains. While it is clear that 
a more elaborate theory is needed to successfully reproduce the 
scattering properties of the grains, it is noteworthy that both the 
compact spheres and the DHS methods predict radial polariza¬ 
tion features (i.e., negatives in the 2^ image) for several types of 
grain distributions, and an increment in signal in the images 
with multiple scattering. 


4. Conclusions 

We investigated the conditions for which non azimuthal polar¬ 
izations can be produced at significant levels in dusty circum- 
stellar environments, and in particular in young protoplanetary 
disks. By doing so we investigated as well the validity of the 
assumption currently made when reducing PDI data using the 
radial Stokes parameters, namely that the linear polarization in 
disks is strictly azimuthal and any non-azimuthal polarization is 
produced by noise and/or instrumental effects. 

Protoplanetary disks are usually optically thick at NIR wave¬ 
lengths and multiple light scattering alone, or in other words the 
scattering of polarized light, can readily produce detectable sig¬ 
natures in the images. These signatures increase with inclina¬ 
tion angle and optical thickness, and this trend holds for compact 
and porous grains. In the disk models discussed above the signal 
in U(p can reach up to ~ 4.5% of the Pj (unconvolved, peak-to- 
peak) for a disk inclined by 40°. This value goes up to 50% for an 
inclination of 70°for the most massive disk shown in Fig.[^ Fur¬ 
thermore, dust properties alone may also lead to violation of the 
strictly azimuthal linear polarization assumption. Astronomical 
observations, laboratory measurements, and models indicate that 
it is quite common in nature to find dust chemical compositions, 
shapes and size distributions that will produce azimuthally cen- 
trosymmetric polarization for most scattering angles but that will 
rotate the polarization plane by 90° for large scattering angles. 
This will produce radial polarization resulting in negative signal 
in the 2^ images. This is often the case in comets and similari¬ 
ties are found, at least in terms of chemical composition (spec¬ 
troscopic signatures), between cometary dust and protoplanetary 
disk dust. 

We have demonstrated that genuine astrophysical signal in 
the U(p images is expected for a number of very plausible config¬ 
urations in protoplanetary disks and other circumstellar environ¬ 
ments. A detailed analysis of these scenarios and their signatures 
in the 2^ and U(p images is beyond the scope of this letter but we 
stress that caution is in order when assuming that linear polar¬ 
ization is purely azimuthally centrosymmetric in protoplanetary 
disks, in particular for disks with inclinations i > 40°. At lower 
inclinations our models show that the images usually contain 
low signals (at least for the disk models tested here), and it is 
probably safe to use the radial formalism. If significant signal is 
found in the U(p images, the formalism may remain relevant for 
disk detection experiments, but for quantitative studies (surface 
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Fig. 2. From top to bottom: synthetic images at l.ljdm for three models (/ = 40°) with outer disk masses of M^ust 2 = 5 x [10 10 10 Mq. 

From left to right, unconvolved images, and convolved t/^, Pj. All models were computed using the DHS method and share the same grain 
populations and porosity (see text). For each row, the images are scaled to the P/ image with the maximum value of Pj being 100. Neither noise 
nor instrumental polarization are included in these simulations, highlighting the relative changes in the and images with multiple scattering. 


brightness, colors, polarization levels) we recommend to analyze 
as well the “traditional” Q and U images. On the positive side, 
signal in the images once all instrumental and data-reduction 
artifacts have been corrected for implies that there is more infor¬ 
mation to be extracted regarding the dust population and particle 
density. 
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